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- - - We discuss three applications of a gauge theory of gravity to rotating as- 

\ trophysical systems. The theory employs gauge fields in a flat Minkowski 

Q"^ ' background spacetime to describe gravitational interactions. The iron fluo- 

rescence line observed in AGN is discussed, assuming that the line originates 
from matter in an accretion disk around a Kerr (rotating) black hole. Gauge- 
theory gravity, expressed in the language of Geometric Algebra, allows very 
• efficient numerical calculation of photon paths. From these paths we are 

able to infer the line shape of the iron line. Comparison with observational 
data allows us to constrain the black hole parameters, and, for the first time, 
J> ' infer an emissivity profile for the accretion disk. The topological constraints 

^ . imposed by gauge-theory gravity are exploited to investigate the nature of 

the Kerr singularity. This reveals a simple physical picture of a ring of mat- 
' ter moving at the speed of light which surrounds a sheet of pure isotropic 

^ ■ tension. Implications for the end-points of collapse processes are discussed. 

Q>^ ! Finally we consider rigidly-rotating cosmic strings. It is shown that a so- 

I lution in the literature has an unphysical stress-energy tensor on the axis. 

Q^' Well defined solutions are presented for an ideal two-dimensional fluid. The 
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exterior vacuum solution admits closed timelike curves and exerts a confining 



\ force. 

a 

1 Introduction 

X 

■ The problem of formulating gravitational theory clS db £lcIU. ge theory has been consid- 

ered by several authorsS'i. In the previous Erice lecturesH, some of the present authors 
(with Stephen Gull) presented a gauge theory of gravity which employed a pair of 
gauge fields defined over a fiat (structureless) Minkowski spacetime (see Lasenby et 
al. u for a complete treatment). This theory provides a radically different picture of 
gravitational interactions from that of general relativity. Despite this, the two theo- 
ries agree in their predictions over a wide range of phenomena. Important differences 
only start to arise over global issues such as the role of topology and horizons, and 
the interface with quantum theory. 

In this lecture we consider the application of this theory to three astrophysical 
situations involving rotating matter. The first application is to the iron fluorescence 
line from the accretion disk around a black hole. X-ray observations of MCG-6-30-15 
show that the iron lines for this Seyfert-1 galaxy are broad and skew i'i. Fits to 
the line profile suggest that the lines originate from fluorescence of matter from the 
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surface of an accretion disk in the strong gravity region around a rotating black hole. 
Modelling the line profile requires the integration of photon trajectories in the region 
of spacetime outside the horizon. Since we are only concerned with properties outside 
the horizon, the predictions of gauge-theory gravity and general relativity coincide 
here, although the gauge theory approach provides much improved machinery for 
performing these integrations. 

The second application is a study of the nature of the singularity at the centre 
of a Kerr black hole It This application fully exploits the fact that in gauge-theory 
gravity, gravitational interaction is mediated by gauge fields defined over a flat back- 
ground spacetime with trivial topology. By integrating the stress-energy tensor over 
the singular regions we reveal a surprising, but physically simple, structure to the sin- 
gularity. These predictions are quite different from the (maximally extended) solution 
favoured by general relativity. 

Our final application is a brief discussion of rigidly-rotating string solutions!. We 
restrict attention to solutions where the direction along the string axis drops out of 
the dynamics entirely, so that we effectively model gravity in (2+l)-dimensions. The 
solution of Jensen and Soleng0 describing a finite width rotating string falls into this 
class of solutions. However, we show that the stress-energy tensor derived from their 
solution is unphysical since it is ill-defined on the string axis. This problem is easily 
overcome, and we close by presenting a set of analytic solutions for rigidly-rotating 
cosmic strings. 

We have found that the Geometric Algebra of spacetime — the Spacetime Algebra 
(STA)ini — is the optimal language in which to express gauge-theory gravity. Employ- 
ing the STA not only simplifies much of the mathematics, but it often brings the un- 
derlying physics to the fore. We begin with a brief introduction to Geometric Algebra, 
the STA and gauge-theory gravity. We employ natural units {G = c = h = to = 1) 
throughout this lecture, except when expressing numerical results. 

2 Geometric Algebra 

This brief introduction to Geometric (or Clifford) Algebra is intended to establish our 
notation and conventions. More complete introductions may be found in Lasenby et 
alM and Hestenes0. The basic idea is to extend the algebra of scalars to an algebra 
of vectors. We do this by introducing an associative (Clifford) product over a graded 
linear space. We identify scalars with the grade elements of this space, and vectors 
with the grade 1 elements. Under this product scalars commute with all elements, and 
vectors square to give scalars. If a and b are two vectors, then we write the Clifford 
product as the juxtaposition ab. This product decomposes into a symmetric and 
an antisymmetric part, which define the inner and outer products between vectors, 
denoted by a dot and a wedge respectively: 

a-b = ^{ab + ba) 
aAb = j{ab — ba). 
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It is simple to show that a-b is a scalar, but aAb is neither a scalar nor a vector. It 
defines a new geometric element called a bivector (grade 2). This may be regarded 
as a directed plane segment, which specifies the plane containing a and b. Note that 
if a and b are parallel, then ab = ba, whilst ab = —ba for a and b perpendicular. This 
process may be repeated to generate higher grade elements, and hence a basis for the 
linear space. 

2.1 The Spacetime Algebra (STA) 

The Spacetime Algebra is the geometric algebra of spacetime. This is familiar to 
physicists in the guise of the algebra generated from the Dirac 7-matrices. The STA 
is generated by four orthogonal vectors {7/^}, /i = . . . 3, satisfying 

It^-lu = + Ivlt,) = Vfiu = diag(+ ). (2.2) 

A full basis for the STA is provided by the set 

1 {7/.} Wk,icrk} {i'j^} i 

1 scalar 4 vectors 6 bivectors 4 trivectors 1 pseudoscalar (2-3) 
grade grade 1 grade 2 grade 3 grade 4 

where ak = 'Jklo, k = 1 . . .3, and i = 70717273 = o'ia2cr3. The pseudoscalar i squares 
to —1 and anticommutes with all odd-grade elements. The {0"^} generate the geo- 
metric algebra of Euclidean 3-space, and are isomorphic to the Pauli matrices. They 
represent a frame of 'relative vectors' ('relative' to the timelike vector 70 employed in 
their definition). The {ck} are bivectors in four-dimensional spacetime, but 3- vectors 
in the relative 3-space orthogonal to 70. We will often denote relative vectors in bold 
typeface (the {cxfc} being the exception). 

An arbitrary real superposition of the basis elements ( |2.3| ) is called a 'multivector', 
and these inherit the associative Clifford product of the {jfi} generators. For a grade- 
r multivector A^ and a grade-s multivector Bg we define the inner and outer products 
via 

Ar-Bg = {ArBs)lr-s\, ArABg = {ArBs)r+s, (2.4) 

where {M)r denotes the grade-r part of M. We shall also make use of the commutator 
product, 

AxB = 1{AB - BA). (2.5) 

The operation of reversion, denoted by a tilde, is defined by 

(AB) = BA (2.6) 

and the rule that vectors are unchanged under reversion. We adopt the convention 
that in the absence of brackets, inner, outer and commutator products take precedence 
over Clifford products. 
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Vectors are usually denoted in lower case Latin, a, or Greek for basis frame vectors. 
Introducing coordinates {x'^{x)} gives rise to a (coordinate) frame of vectors {e^} 
where = dfj_x. The reciprocal frame, denoted by {e'^}, satisfies e^-e'' = S^. The 
vector derivative V(= d^) is then defined by 

V = e'^a^ (2.7) 

where 5^ = d/dx^. 

Linear functions mapping vectors to vectors are usually denoted with an underbar, 
/(a) (where a is the vector argument), with the adjoint denoted with an overbar, /(a). 
Linear functions extend to act on multivectors via the rule 

liaAbA ■■■Ac) = /(a)A/(6)A • • • A/(c), (2.8) 

which defines a grade-preserving linear operation. In the STA, tensor objects are rep- 
resented by linear functions, and all manipulations can be carried out in a coordinate- 
free manner. 

All Lorentz boosts or spatial rotations are performed with rotors. These are even- 
grade elements R, satisfying RR — 1. Any element of the algebra, M, transforms 
as 

M ^ RMR. (2.9) 

A general rotor may be written as i? = exp(i?/2) where 5 is a bivector in the plane 
of rotation. 

2.2 Gauge-Theory Gravity 

Physical equations, when written in the STA, always take the form 

A{x) = B{x), (2.10) 

where A(x) and B{x) are multivector fields, and x is the four-dimensional position 
vector in the (background) Minkowski spacetime. We demand that the physical 
content of the field equations be invariant under arbitrary local displacements of the 
fields in the background spacetime, 

A{x)^A{x'), x'^f{x), (2.11) 

with f{x) a non-singular function of x. We further demand that the physical content 
of the field equations be invariant under an arbitrary local rotation 

A{x)^RA{x)R, (2.12) 

with R a non-singular rotor-valued function of x. These demands are clearly equiv- 
alent to requiring covariance (form-invariance under the above transformations) of 
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the field equations. These requirements are automatically satisfied for non-derivative 
relations, but to ensure covariance in the presence of derivatives we must gauge the 
derivative in the background spacetime. The gauge fields must transform suitably 
under (local) displacements and rotations, to ensure covariance of the field equations. 
This leads to the introduction of two gauge fields: h{a) and Q{a). The first of these, 
h{a), is a position-dependent linear function mapping the vector argument a to vec- 
tors. The position dependence is usually left implicit. Its gauge-theoretic purpose 
is to ensure covariance of the equations under arbitrary local displacements of the 
matter fields in the background spacetime iS. The second gauge field, f2(a), is a 
position-dependent linear function which maps the vector a to bivectors. Its intro- 
duction ensures covariance of the equations under local rotations of vector and tensor 
fields, at a point, in the background spacetime. 

Once this gauging has been carried out, and a suitable Lagrangian for the matter 
fields and gauge fields has been constructed, we find that gravity has been intro- 
duced. Despite this, we are still parameterising spacetime points by vectors in a fiat 
background Minkowski spacetime. The covariance of the field equations ensures that 
the particular parameterisation we choose has no physical significance. The feature 
that is particularly relevant to this lecture is that we still have all the features of 
the fiatspace STA at our disposal. A particular choice of parameterisation is called 
a gauge. Under gauge transformations, the physical fields and the gauge fields will 
change, but this does not alter physical predictions if we demand that such predictions 
be extracted in a gauge-invariant manner. 

The covariant Riemann tensor TZ[aAb) is a linear function mapping bivectors to 
bivectors. It is defined via the field strength of the Q{a) gauge field: 

nh-\aAb) = a-Vn{b) - 6-Vfi(a) + n{a)xn{b). (2.13) 

The Ricci tensor, Ricci scalar and Einstein tensor are formed from contractions of 
the Riemann tensor: 

Ricci Tensor: 7^(a) = 7^■7^(7^Aa) (2.14) 
Ricci Scalar: 7^ = 7^■7^(7^) (2.15) 
Einstein Tensor: Q{a) = lZ{a) — \alZ. (2-16) 

The Einstein equation may then be written as 

g{a) = KT{a), (2.17) 

where T{a) is the covariant, matter stress-energy tensor. The remaining field equation 
gives the fi-function in terms of the /i-function, and the spin of the matter field 
However, this will not be required for this lecture. 

Some comments on gauge-theory gravity are now in order. Firstly, we note that 
the theory is formally similar in its equations (hence local behaviour) to the Einstein- 
Cartan-Kibble-Sciama spin-torsion theory i, but it restricts the Lagrangian type and 
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the torsion type [TZ"^ terms in the gravitational Lagrangian, or torsion that is not 
trivector type, leads to minimally coupled Lagraneians giving non-minimally coupled 
equations for quantum fields with non-zero spina). As an interesting aside, we note 
that self-consistent homogeneous cosmologies, based on a classical Dirac field, require 
that k = (the universe is spatially flat) EJ. 

If we restrict attention to situations where the gravitating matter has no spin, 
then there are still differences between general relativity and the theory presented 
here. These differences arise when time reversal effects are important {e.g. horizons), 
when quantum effects are important, and when topological issues are addressed. For 
example, there is no analogue of the Kruskal extension of the Schwarzschild solution 
in our theory. These differences arise from the first-order derivative nature of the 
theory, and its origin in a fiat background spacetimei. 

Even in those cases where the gauge-theory predictions are completely in ac- 
cord with general relativity (all present experimental tests), we believe that our ap- 
proach offers real computational advantages over conventional methods. The 'Intrin- 
sic method' described in Lasenby et a/. 00 is a good example of the power of the 
gauge-theory approach. This method allows the field equations to be solved in vari- 
ables which are covariant under displacement gauge transformations. The first-order 
'rotor' approach to calculating photon trajectories, discussed in the next section, is 
another such example. 

3 The Iron Fluorescence Line 

The X-ray emission from AGN is believed to originate on an accretion disk around 
a black hole. In particular, if the disk material absorbs continuum radiation with 
energy > 7.2 keV, then a fluorescent iron line at 6.4 keV may result (the probability 
for this absorption is high, ~ 0.34 per incident photon). Such lines were observed 
by Pounds et al.^ and Matsuoka et a/. in Seyfert-1 galaxies. Recent observations 
of MCG-6-30-15 {z = 0.008) have shown that this line is both broad and skewl'l. 
Figure |I| shows the line profile from Iwasawa et al. i, which is averaged over the 
1.7 X 10^ s observation period, and normalised to a power-law model (which included 
corrections for cold absorption). The broad iron K emission line lies around 6keV. 
Recent work on the variability in the line profile during the observation has shown that 
the line shape varied with position on the light curve (see Figure ^ for the observed 
light curve, reproduced from Iwasawa et a/.S) and that at the minimum emission, the 
lineshape broadened further. In particular, the lineshape extended further to the red 
side and the blue wing disappeared. The line flux at minimum emission is shown in 
Figure |^, which should be compared to the average line flux over the entire observation 
(Figure 0). The redshift factor at the tail of the red wing extends to around ~ 0.5, 
showing that we are seeing the effects of very strong gravity at the epoch of minimum 
emission. If this redshift were due to climbing out of a Schwarzschild (non-rotating) 
black hole, then the emission would have to occur from r ~ 2.5GM/c'^, where M is 
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Figure 1: The ratio of data and model for the averaged 0.4-10 keV spectrum of MCG-6-30-15. The 
data are obtained by integrating over the entire observation (~ 1.7 x 10^ s). The model is a single 
power-law with photon index 1.96, modified by cold absorption, fitted to the data excluding the 
0.7-2.5 keV and the 4.5-7.2 keV bands. There is a clear absorption feature around IkeV due to a 
warm absorber, and a broad iron K emission line around 6keV. Reproduced with permission from 
Iwasawa et al. 

the mass of the black hole. However, the minimum radius stable circular orbit in a 
Schwarzschild black hole is at 6GM/ (? . For a Kerr (rotating) black hole this minimum 
radius goes down to GMjc? for a corotating orbit. The most likely conclusion is that 
the black hole is rapidly rotating. 

We shall assume that the variability in line profile is due to flaring and that at 
minimum emission, we are seeing only the effects of a uniform accretion disk. Previous 
authors llilli have calculated the predicted lineshape for a maximal Kerr black hole, 
but in order to fit the lineshape properly we must predict the lineshape for arbitrary 
angular momentum, inclination angle (angle between the line of sight of the observer 
and the axis of rotation) and accretion disk parameters. This problem was addressed 
by a collaboration including two of the present authors 0. 

3.1 Predicting the lineshapes 

In order to predict the lineshape, we require the redshift and point of intersection 
with the accretion disk, for all those null geodesies passing through the observation 
point and the accretion disk (in the past). Gauge-theory gravity is particularly useful 
here, since we can employ a computationally efficient 'rotor' approach to the prob- 
lem. This approach arises naturally in several diverse settings, including the motion 
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Figure 2: The 0.5-10 keV light curve from MCG-6-30-15. The epoch of the start of the hght curve 
is 1994 July 23 05:05:25. Each data bin is averaged over 128s. Reproduced with permission from 

Iwasawa et al. 
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Figure 3: The observed iron-line flux from MCG-6-30-15 at minimum emission. 
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of charged particles in electromagnetic fields ta, and the motion of particles in grav- 
itational fields (including torsion effects) 0. The rotor approach is useful not only 
because of the computational efficiency of the resulting first-order equations, but also 
because of their numerical stability. We have found that these first-order techniques 
are generally faster and more accurate than direct integration of the (second-order) 
geodesic equations. 

We begin by parameterising the photon 4-momentum with the aid of two rotors, 

where a and f3 are scalar functions of the affine parameter A along the null geodesic. 
We then form the rotor R = R1R2, which directly controls the direction of the photon 
4-momentum p via 

p = $i?(7o + 71)^, (3.2) 

where $ is another scalar function of A, which equals the energy of the photon relative 
to an observer with covariant 4- velocity 70. Note that p is guaranteed to be null since 
7o + 7i is null. 

The basic dynamical equations are@ 

p = —Q{x)-p (3.3) 
x = h{p), (3.4) 

where x is the spacetime position vector of the photon, and overdots denote differen- 
tiation with respect to A. For the /i-function we use the form appropriate to the Kerr 
black hole in Boyer-Lindquist form: 



= — ^TTPT^t ^</" ^(^r) = er 

pA^'^ rp p 

— ar"^ ?AV?6 r r 

h[e^) = . et + -60, h{eg) = -eg, 

pA^/^ p p 



(3.5) 



where 



p = + cos^^, A = r^-2Mr + a^, (3.6) 



a is the black hole angular momentum, and M is its mass. The vectors appearing 
in ( pj.5D are the polar frame vectors associated with the polar coordinate system 
{t,r,^,0}: 

t = x-jo cos6 = x-'j^/r 

(3.7) 



r = (xA7o)^ tan0 = (x-7^)/(x-7^). 
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The /i-function given by ( |3.5| ) is singular where A = 0. This therefore fails to define a 
global solution. A global solution can be obtained by a (singular) gauge transforma- 
tion. The resulting solution would allow discussion of properties inside the horizon 
(see Section although the above form is adequate to describe the spacetime exterior 
to the horizon. The Riemann tensor associated with (13.51) takes the neat form 



-M 

^(^) = 777 : 7m + Se.e^^e.et), (3.8) 

which is clearly non-singular over its domain of validity. 

The model described here cannot discriminate a and M separately. Instead the 
relevant black hole parameter is a^, = a/M. For an extreme Kerr black hole a^, = 
1, whilst the most extreme stable system of hole and accretion disk probably has 
a=K=0.998lll. The remaining parameters to which the line profile are sensitive are the 
inclination angle i (this double usage of the symbol i should not cause any confusion), 
and the radial emissivity profile e(r). 

3.2 Numerical results 

The equations for p and x yield seven first-order differential equations in A for the 
quantities $,a; and (3 (giving the photon 4-momentum), and the coordinates of the 
photon t, r, 6 and (f). These equations are not only easy to solve numerically, but 
automatically conserve energy, angular momentum and Carter constant (see, for ex- 
ample. Section 33.5 of Misner et alx^ for a discussion of the (fourth) Carter constant). 
For each photon path from the accretion disk to the observer, the redshift may be 
calculated using 

l + z = = , (3.9) 

^obs ^obs 'Pobs 

where Vem and fobs are the covariant 4-velocities of the emitting gas and a distant 
observer respectively, and Uem and Uohs are emitted and observed frequencies. We 
assume that the matter responsible for the fluorescent line lies on the (geometrically 
thin) accretion disk and has velocity f em given by the velocity of a circular equatorial 
geodesic. Figure § shows the predicted frequency contrast of the iron line in the 
images of the disk, as seen by the distant observer. Light bending due to the strong 
gravity near the black hole is clearly visible in the images at large inclination angles. 

The fluorescent line is emitted with intensity /em(i^em) in the local rest frame of the 
disk, and received with intensity /obs(^'obs) by the distant observer. These intensities 
are related by the invariant along the photon path 

-^obs ( '^obs ) -^om ( ^em ) /" Q i n ^ 

obs em 

We may then integrate over the solid angle subtended at the observer, and the fre- 
quency bins of the detector to obtain the flux. The emitted intensity lemi^em) will be 
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Figure 4: Images of the accretion disk observed in the iron-hne by a distant observer. The top row 
of images is for a Schwarzschild black hole (a* = 0) . The middle row shows the images for a Kerr 
black hole with a* = 0.5. The bottom row corresponds to an extreme Kerr hole (a* = 0.998). For 
each value of a*, the image is shown for inclination angles 30°, 75° and 90°. The grey-scale indicates 
the variation of redshift in the image. 
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a function of the radius of emission in the accretion disk. Typical assumptions are 
that the disk continuum emissivity follows e(r) oc r~'^ with q in the range 2-3, or that 
e(r) follows the (more realistic) law of Page and Thorne0. We make the assumption 
that the line emissivity follows the latter, with emission starting at the radius of 
marginal stability. In Figure ^ we show nine predicted line profiles for different values 
of a=K and inclination angle. These suggest that for fixed emissivity profile, the overall 
line shape is most sensitive to the inclination angle. For sufficiently large inclination 
angles the spectrum is double peaked. This effect is mainly due to Doppler shifts 
between the receding and approaching parts of the disk0. The two-parameter {a^ 
and i) predicted lineshapes may now be fitted to the observed profile over the period 
of lowest luminosity. The confidence contours in the parameter space are shown 
in Figure ^ along with the best-fit to the observed line fiux. The contours give 
relatively strong evidence for an inclination angle of ~ 25° - 30°, and strongly favour 
an extreme Kerr black hole (a* > 0.94). 

The need for a high value of is so crucial to fitting the observed line profile, 
that the observed line strongly constrains the emissivity profile e(r). This allows us 
to infer an emissivity profile for the first time. We have done this with the other 
parameters a^, and i held fixed at the values 0.998 and 30°. The inferred profile is 
shown in Figure |^. There is some evidence for a power law, with a value of g ~ 3.5, 
although it obviously becomes noisy at low fiux levels. 

4 The Nature of the Kerr Singularity 

The problem we now wish to address concerns the endpoint of rotating collapsing 
matter. In the previous section we assumed that the collapsed system at the centre 
of the AGN was a Kerr black hole. For this reason it is of interest to look at the nature 
of the singularity inside a Kerr black hole, and so determine whether it is consistent 
with what we would expect for the end point of such a collapse. This problem was 
considered in Doran0 and Doran et al. i. 

Gauge-theory gravity allows an unambiguous answer to this problem, provided 
that one accepts the basic premises of our theory. This is because the topological 
constraints implied by the theory ensure that we have well defined surfaces over 
which we may apply integral theorems. The results of this investigation will be a 
gauge invariant description of the nature of the singularity, but so far we have found 
that the problem is most tractable with a specific choice of gauge. We shall employ 
the 'Kerr-Schild' gauge, in which the /i-function takes the form 

h{a) = a + a-U, (4.1) 

where / is a null vector [P = 0). This /i-function is globally valid, unlike the gauge 
employed in the previous section. This is essential to study the nature of the sin- 
gularity, since this lies inside the horizon. A simple example of a solution in the 
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Figure 5: Calculated iron- line profiles as measured by the distant observer. The top row is for 
a* = 0.998, the middle row has a* = 0.5, and the bottom row is for a Schwarzschild black-hole 
(a* = 0). For each value of angular momentum, line profiles are given for inclination angles of 
0°,30° and 89°. 
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Figure 6: Figure 6a shows the contours of probability versus i and a*, calculated assuming a uniform 
prior on each parameter. Figure 6b compares the calculated line flux (dashed line) corresponding 
to the best-fit values i = 29°^1'2, a* = l^o oi with the observed data (marked points). The error 
bounds on the best-fit parameters are 68% confidence intervals obtained by marginalization over the 
other parameter. 
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Figure 7: A log-log plot of the inferred emissivity profile, for a* = 0.998 and inclination angle 
i = 30°. 



Kerr-Schild form is provided by a Schwarzschild black hole, which has 
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\/m7^(7o - er 



(4.2) 



In this gauge, incoming radial photons follow straight lines in a {t,r) plot, and ter- 
minate on the singularity at r = 0. Outgoing radial photons may only escape from 
the black hole if they start outside the horizon (which lies at r = 2M). This solu- 
tion is geodesically incomplete and is not time-reverse symmetric. This forces us to 
adopt the picture of the black hole being the end-point of a collapse process, with the 
formation of the horizon capturing information about the direction of time for which 
the collapse occurredifl 

It is easy to show that for a general Kerr-Schild vacuum solution!. 



/■V/ oc /. 



(4.3) 



We shall only consider matter fields for which this relation is also true (this clearly 
restricts the matter fields that we may describe, but does include the Reissner- 
Nordstrom and Vaidya 'shining star' solutions). It follows that we may write 



(4.4) 



where (p is an arbitrary scalar function of position. The Einstein tensor then takes 
the form 



(4.5) 
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where 



uj{a) = VA{a-ll) (4.6) 



and the vector a is not differentiated. Equation ( [4.5|) shows that for this class of fields, 
the Einstein tensor is a total divergence in the background spacetime. For stationary 
fields, it follows that the Einstein tensor is a total 3-divergence. This allows us to 
convert integrals of Q{a) over the singular regions of space to surface integrals over 
well defined 2-surfaces enclosing the singularity. 



For example, the Schwarzschild black hole (^4.2|) gives 



/ d^xGia) = 8nMa-fo-fo, (4.7) 

J r<ro 

where ro is any value > 0, since Q{a) vanishes everywhere except at the origin. It 
follows that the matter stress-energy tensor is given by 

r(a) = MS{x)a-folo, (4.8) 

where x = xA'-fo. This is the stress-energy tensor appropriate to a point source of 
matter (of mass M) following the world line r = 0. This technique is analogous 
to the usual analysis of the singularity in the Coulomb field, due to a point charge. 
Note that the integrals that we have performed are not gauge invariant, but we have 
extracted gauge covariant information in the form of the stress-energy tensor. 

4-1 The Reissner-Nordstrom solution 

We now highlight a result obtained by one of us in Doranlzl. The Reissner-Nordstrom 
solution describes a charged, non-rotating black hole. In the Kerr-Schild gauge, the 
solution may be written in the form 

/i(a) = a + ?7a-e_e_, (4.9) 

where, 



M 

Svrr 



V = — -^ZZ^^ e_=7o-e^, (4.10) 



and q is the charge of the source. Away from the origin, the stress-energy tensor 
evaluates to 

T(a) = -\TaT, (4.11) 

where T = ger7o/(47rr^). This is the expected form for the electromagnetic stress 
tensor due to a point charge q at the origin. To study the behaviour in the singular 
region, we return to equation ( [4.5| ) to obtain 

2 

/ d^x T (a) = Ma-7o7o + 777— (« - 4a-7o7o) • (4. 12) 

Jr<ro IAtXTq 
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The first term on the right-hand side is the same as in the Schwarzschild case, whilst 
the latter is the trace-free electromagnetic contribution. Concentrating on the 70- 
frame energy component, we find that 



2 

/ ci=^x7o-r(7o) = M--^. (4.13) 

Jr<ro ovrro 



Something remarkable has happened here — due to the gravitational fields,, the elec- 
tromagnetic contribution to the energy is now negative and vanishes as we extend the 
integral over all of space {tq — >■ 00). This is in stark contrast to the standard picture 
from classical electromagnetism, where the self-energy of the point charge diverges. 
Inclusion of the gravitational fields has removed this divergence, ensuring that the 
total electromagnetic self-energy vanishes. The manner in which this regularisation 
is achieved is discussed in Doranl. 

4-2 The Kerr solution 

The Kerr solution describes a rotating, uncharged black hole. A remarkable complex 
harmonic structure underlying this solution was found by Schiffer et al. 11. We define 
'complex' numbers 7 and uj via 

7 = a-Fi/3, = 7"\ (4.14) 

where a and (3 are scalars. Not that the 'z' appearing in equation ( [4.14| ) is the 
spacetime pseudoscalar. This element is the generator of duality transformations. 
For example, the STA statement of the self-duality of the Weyl tensor is 

W{iB)=iW{B), (4.15) 

where B is an arbitrary bivector. We obtain an axisymmetric solution of Kerr-Schild 
form if we can solve the two equations 

V^7 = 0, {Vuf = 1, (4.16) 

where V = 70AV is the derivative operator in the space orthogonal to 70. Any 
solution of these equations generates a Kerr-Schild type solution of the form 
with / given in terms of 7 and 70 i. 

As a simple example, the Schwarzschild solution is obtained by setting u = r. We 
can obtain the Kerr solution by a 'complex translation' of the Schwarzschild solution: 

u;=(x^ + y^ + (z-tL)Y^ (4.17) 

where L is a scalar constant, and {x,y,z} are Cartesian coordinates. The Riemann 
tensor for this solution evaluates to 

M 

n{B) = -^{B + Sa^Ba^), (4.18) 
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with the unit bivector cr^ given by 

X — Lia^i 

= -. 4.19 

The Riemann tensor is only singular where o; = which occurs on the ring p = 
L, z = {p = (x^ + 1/^)^/^). For this reason, it has been widely believed that the Kerr 
singularity is a ring only. 

We can analyse the nature of the Kerr singularity in a similar manner to the 
Schwarzschild and Reissner- Nordstrom cases treated earlier. We begin by integrating 
over a spatial region which fully encloses the central disk. We find that 

J d^xT{a) = Ma-7o7o, (4.20) 

where M is the mass of the hole (this constant appears when relating u to I), and the 
integral is taken over any region enclosing the central disk. This is the same result 
as in the Schwarzschild case. We can also integrate the (orbital) angular momentum 
tensor xAT{a) {x being the spacetime position vector) over the region enclosing the 
disk to obtain 

J d^xxAT{a) = ML[-a-7oicr3 + |(aA7o) xicrg]. (4.21) 

This clearly identifies ML as the total angular momentum in the fields, as expected 
from their long-range behaviour. 

To examine the matter distribution for p < L, we integrate the Einstein tensor 
over cylindrical 3- volumes normal to the disk. The calculations are lengthy and great 
care must be taken over the choice of branch for the complex square roots. Details 
are given in Doran et a/, i, where it is shown that for p < L, 

6^(73) = 0. 

The vectors Cp and are unit spacelike basis vectors in the cylindrical polar coordinate 
system. This form for the Einstein tensor clearly shows that matter is not located 
solely on the ring at p = L, but also over a disk in the plane z = 0, which has the 
ring as its boundary. We see immediately that T(a) is symmetric, showing that there 
are no hidden sources of torsion in the disk. This contribution to the Einstein tensor 
describes a rigidly-rotating, massless disk of pure isotropic tension in the plane of 
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the disk. The tension is given by M/[47rL(L^ — p^Y^'^]. The angular velocity is 1/L 
so that the edge of the disk follows a lightlike trajectory. Remarkably, this tension 
field has a simple non-gravitational explanation. The special-relativistic equations 
governing a massless, rigidly-rotating membrane (with a ring of particles attached to 
the edge) reproduce exactly the functional form with p just found for this tension. 
The fact that the disk has vanishing energy density but generates a tension means 
that it violates the weak energy condition. 

The integral of the Ricci scalar over the interior of the disk yields SttM, which 
is equal to the value deduced from integrals enclosing the entire singular region. It 
follows that any matter in the ring at p — L makes no contribution to the Ricci scalar, 
and hence that the contribution to the stress-energy tensor from the ring singularity 
must have vanishing trace. Furthermore, the disk of pure isotropic tension can make 
no contribution to the angular momentum of the fields, so the angular momentum 
must come solely from the ring singularity. From these considerations, we may deduce 
that the matter in the ring follows a light like trajectory. These conclusions are gauge 
invariant, since they are inferred from the eigenvalue structure of covariant tensors. 

We see that within the framework of gauge-theory gravity, the Kerr singularity is 
composed of a ring of matter, moving at the speed of light, which surrounds a disk 
of pure isotropic tension. The tension in this disk has precisely the form expected on 
the basis of special-relativistic arguments. The rotating ring of matter is a perfectly 
satisfactory endpoint for matter collapsing with angular momentum — the proper ra- 
dius coincides with the minimum size allowed by special relativity, for an object with 
angular momentum ML. However, the presence of the disk of tension is problematic 
— no baryonic matter can have a tension but vanishing energy density. If baryonic 
matter cannot form this disk, then what is the status of the Kerr solution as the end- 
point of the collapse process? The answer to this question must await the discussion 
of realistic collapse processes within the framework of gauge-theory gravity. 

5 Rigidly-Rotating Cosmic Strings 

As a final topic, we shall turn to a situation with cylindrical symmetry. We shall 
restrict attention to string solutions in which the direction along the string axis plays 
no part in the dynamics of the string. Imposing this restriction means that the 
solutions which include pressure will violate the boost invariance, which is usually 
demanded of all cosmic string solutions. However, these solutions may still be of use 
for rotating strings in (3+l)-dimensions, where it is not clear that one can impose 
boost invariance. 

We adopt a cylindrical polar coordinate system {t, p, 0, z}: 



The vectors {et, Cp, e<^, e^} comprise the associated coordinate frame, with et = 70 



p = ^[-{xAasf] 



t = x-7o 
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and Cz = 73. The reciprocal frame vectors are denoted as {e*, e'', e''^, e^}. We shall 
consider solutions described by an /i-function of the form 

h{e') = he' + p/^e^ /^(e^ = 

We require that the /i-function (and the ^-function) be well defined on the string 
axis {z = 0). This requires that f 21 phi and /12 all vanish smoothly on the axis. 
These requirements replace the notion of 'elementary flatness' employed in the gen- 
eral relativity literature El. This is an area where gauge-theory gravity offers clear 
advantages over general relativity — since we deal solely with linear functions defined 
over a (flat) background spacetime, there is never any doubt about the conditions that 
these functions should satisfy. 



5.1 The solution of Jensen and Soleng 

The first published solution describing the interior of a finite width rotating string 
was that of Jensen and SolengS. Their solution may be generated from an /i-function 
of the form (15.21) with 



1 M 

hi = — coshw fi = cosh-u sinhw 

A A 

ho = — sinhw fo = sinhw sinhw, 

A A ' 



(5.3) 



where 



and 



A = ^sm{VXp), (5.4) 



M = 2a (^{p - ps) cos{y/Xp) - sm{^/Xp) + psj . (5.5) 

Here A is a positive constant, a is a constant with a < 1, and ps is the radius of 
the string. The parameter u appearing in ( p.3| ) is arbitrary up to the constraint that 
u = on the axis of the string (so that the /i-function is well defined there). Analysing 
the solution in the gauge in which u = everywhere, we find that 

Q{et) = -a2et + aX(f), (5.6) 
where a2 is a function whose explicit from we do not require. The vector (p is given 

by 

4) = e^/p = - sin07i + cos072. (5.7) 
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If we now consider an observer with covariant velocity et passing through the axis 
of the string, it is clear that the 3-momentum density he measures on the axis is 
ill-defined. 

The conclusion is that the solution of Jensen and Soleng does not define a physi- 
cally acceptable matter distribution. This is surprising, since the solution does satisfy 
the criteria of elementary flatness. This point illustrates a further advantage of the 
gauge theory approach over general relativity; the gauge theory focuses attention on 
the physically relevant quantities, such as the eigenvalues of the stress-energy ten- 
sor. In such an approach, it quickly becomes apparent if a solution has unphysical 
properties. 



5.2 Rigidly-rotating strings 

It is not difficult to find rotating string solutions with a physically acceptable matter 
distribution. The simplest model is that of a two dimensional ideal fiuid, with stress- 
energy tensor 

T(et) = net 



T{e^) = {fi- 2p)e^. 

The covariant 4-velocity of the fiuid is et, the energy density is fi{p) and the isotropic 
pressure in the ia^ plane is p{p). The coefficient of Q{ez) is restricted to /i — 2p by the 
Einstein equations, and the assumed form for the /i-function (|5.2|) . This stress-energy 
tensor is well defined on the axis, provided that the pressure and energy density are 
finite there. 

A rigidly-rotating (shear-free) solution is given by 
hi = — — /: 



sinAp A + cosAp 



2 ,^ (5.9) 



X{l + A) sinAp' 

where A is an arbitrary positive constant, and the constant A satisfies A < —1. It is 
simple to show that this linear function is well defined on the axis. The pressure and 
energy density evaluate to 

87rp = K^-GT (5.10) 
87ifi = 3K^ + X"^, (5.11) 

where the functions G, K and T are given by 

smAp [A + cosApj^ A + cosAp 
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with B a further constant. The functions G, K and T arise naturally in the rotation- 
gauge field. The boundary of the string occurs where p = 0, and this must be reached 
before p > n/X. 

This interior solution matches onto the exterior vacuum solution given by 

/i = -(1 + A)ia/BY/'[{p + p,f - a\^l^ (5.13) 

= ia/BYf^X^l^P±pl^p^ (5.14) 
(P + Po) 

/2 = 77^^(/i' - 1) (5-15) 
/i(P + Po) 

h2 = 0. (5.16) 

The constants po and a must be determined by the matching conditions at the string 
boundary. This class of vacuum solution does not appear to correspond to anything 
given previously in the literature. There is a confining force in the vacuum meaning 
that no particle can escape from the string, regardless of the initial velocity that it is 
given. 

The line element associated with this external solution is 



a{l + A) 



This class of rigidly-rotating strings is of particular interest because these solutions 
always admit closed timelike curves at some distance from the string. This follows 
from the fact that for large p, fi varies as 1/p whereas /2 tends to a constant value. 
Beyond the point where the magnitude of /2 exceeds that of /i, a closed circular path 
around the string becomes timelike. The long-range properties of these solutions make 
them ultimately unphysical, but there is no reason to suppose that the solutions will 
not be relevant near a string of finite length. 
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